DTP/97/68 
UCSD/PTH 97-16 
hep-ph/9707496| 
July 1997 



Analytic Calculation of Two Loop Corrections to Heavy Quark 
Pair Production Vertices Induced by Light Quarks 



A.H. Hoang a , T. Teubner' 



a Department of Physics, University of California, San Diego, 
La Jolla, CA 92093-0319, USA 

b Department of Physics, University of Durham, 
Durham, DH1 3LE, UK 



Abstract 

Non-singlet 0(a 2 s ) corrections from light quarks are calculated analytically for the massive 
quark-antiquark pair production rates through the vector, axial-vector, scalar and pseu- 
doscalar currents for all energies above the threshold. The results are presented in terms of 
the moments of the absorptive part of the vacuum polarization function of the light quarks, 
which allows for an immediate determination of 0{a 2 s ) corrections coming from other types 
of light particles. The kinematic regime close to the threshold is examined, and a compari- 
son to results based on non-relativistic considerations is carried out. 
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1 Introduction 



In the past years considerable amount of work has been invested to obtain higher order QCD radiative 
corrections to quark pair production vertices. As far as the quark pair production rates through the 
vector and axial-vector currents are concerned this effort was mainly motivated by the impressive 
amount of data on hadronic Z decay events obtained at the LEP experiments. Whereas the 0(a s ) 
corrections have been well known since quite a long time for all energy and mass values @], the 
complete analytic form of the 0{a 2 ) corrections is known as a high energy expansion in M 2 /q 2 
where yfq 2 denotes the cm. energy and M the mass of the produced quarks. For high energies, 
even the leading leading O(a^) @] and some mass corrections have been calculated ||. The quark 
pair production rates through the scalar and pseudoscalar currents have also been investigated in 
view of the expected studies of hadronic Higgs decays at forthcoming collider experiments. Here, the 
0(a s ) corrections are well established for all ratios M 2 /q 2 || [?]], while the second order corrections 
are known as an expansion in M 2 /q 2 [H]. Thus, to 0(a 2 ) the present calculations are definitely 
a good approximation if the cm. energy is much larger than the mass of the produced fermions, 
which is for example the case for Z or Higgs decays into bottom quarks, but can also be applied 
for intermediate energies if a sufficient number of terms in the high energy expansion is taken into 
account (see also M for a review). However, in face of forthcoming experiments (r-charm-, B-factory, 
NLC), where quark pairs will also be produced in the kinematic regime where the quark masses and 
the cm. energy are of comparable order, i.e. closer to the threshold, the knowledge of the entire 
0(a 2 ) corrections valid for all values of M 2 /q 2 is desirable. Whereas meanwhile there are methods 
to obtain numerical approximations for the 0(a 2 ) corrections in the latter kinematic regime based on 
Pade approximants Eo|, the entire analytic calculation seems to be an impossible task at the present 
stage. However, analytical results are of particular importance because they provide important cross 
checks for approximation methods and serve as a starting point for qualitative considerations. 

A first step towards a complete analytical calculation consists of the determination of those 0(a 2 ) 
corrections to the heavy quark production vertices which are induced by light quarks. These corrections 
are equal to the sum of all possible cuts of at least two massive quark lines of current-current correlator 
diagrams consisting of a massive quark loop, a massless quark loop and two gluon lines. The individual 
contributions can be divided into three classes: (i) the massive quarks are produced directly by 
the external current insertion (primary production), (ii) the massive quarks are produced via gluon 
emission off the massless quarks (secondary production) and (hi) the interference between amplitudes 
of class (i) and (ii). Class (iii) corresponds to cuts of the double-triangle or singlet diagrams. For the 
vector current correlator the class (iii) diagrams vanish due to Furry's theorem. The contributions 
belonging to class (ii) were determined analytically in [jll|, [TJ| for the vector and axial-vector case.Q 

In this work we present analytic expressions for the 0(a 2 ) corrections belonging to class (i) for 
arbitrary values of the quark mass M and cm. energy ^/q 2 above threshold. Analytic formulae for 
primary production through the vector and scalar currents were already presented in []l4j and |lq| , 
respectively. Here we carry out the same program for the axial-vector and pseudoscalar case. In 
addition, we present the results in terms of moments of the absorptive part of the vacuum polarization 
of the light quark-antiquark pair, proceeding along the lines of [16|, which allows for an immediate 
determination of the 0(a 2 ) corrections coming from light colored scalar particles or from the gluonic 

1 A numerical approach for the contributions belonging to class (ii) in the vector and axial-vector case was presented 
in [El. 
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vacuum polarization insertion. Based on the concept of the moments even a certain class of O(af) 
and 0(a A s ) corrections can be determined. Thus, for completeness also the vector and scalar cases are 
presented in terms of the moments. 

To be definite we consider the light quark induced 0(a 2 s ) contributions to the imaginary part of 
the massive quark current-current correlators 

= -iJdxe l ^(0\T3°(x)3?(0)\0) 

= g^q 2 n @ (q 2 )- qtl q u Ut(q 2 ), Q = V,A, (1) 
for the vector and the axial- vector currents 

# = *2,7„*g. tf = *&7,.T6*Q, (2) 

and 

q 2 U e (q 2 ) = iJdxe^ x (0\Tj e (x)j @ (0)\0), Q = S,P, (3) 

for the scalar and pseudoscalar currents 

j s = M0*° Q 9° Q , / = ^°$° 75 *», (4) 

where summation over color degrees of freedom is understood, denotes the bare massive quark 
field and M° its bare mass. In order to maintain vanishing anomalous dimension for all four currents 
the scalar and pseudoscalar currents are defined including the mass parameter. The longitudinal 
part of the vector current correlator II^ is equal to H v due to current conservation, whereas II^ is 
related to the pseudoscalar current correlator via the axial Ward identity. The optical theorem 
relates the renormalization group invariant imaginary parts of II , B = V, A, S, P, to decay rates and 
cross-sections. The contribution of Imll^'" 4 to the total hadronic Z decay rate reads 

AIT = ^^(v 2 Q R v (M 2 ) + a 2 Q R A (M 2 )), (5) 
V Q = 2I 3 - 4 <7Q sin 2 ^, 



where 



R Ay {q 2 ) = 12 7rImn A ' y (^+ze), (6) 



qQ is the massive quark electric charge, @w the Weinberg angle and ig the massive quark weak isospin. 
Moreover, R v is equal to the cross-section of QQ production in e + e~ annihilation into a single photon 
normalized to the point cross-section. The contribution of the quantities R V ' A to Z-mediated QQ 
production is straightforward. The contributions of ImII s '' P to the hadronic decay widths of a scalar 
or pseudoscalar Higgs boson with mass ms and mp, respectively, read 

A -phad Gf M 2 f s p 2 

where 

R S > P (q 2 ) = ™l m ll s > p ( q 2 + ie) (8) 
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and M denotes the heavy quark pole mass. The coefficient a€ P is the relative weight of the cor- 
responding partial width with respect to the standard model Higgs decay rate (a# = 1). Thus 
in the framework of the MSSM, where a two doublet Higgs sector is realized, a^f = sin 2 a/ sin 2 /3, 

a down = cog 2 a j cog 2 ^ a up = cog 2 a / ^ ^ a down = gin 2 a j cog 2 ^ Q up = mt 2 p and a down = ^2 ^ 

where H, h and A denote the heavy and light scalar and the pseudoscalar Higgs fields, respectively, a 
describes the mixing between the two Higgs doublets and (3 is defined via the ratio of their vacuum 
expectation values, tan/3 = V2/V1 (see e.g. (lTj). In the following sections we present the results of our 
calculation in terms of the quantities R® , @ = V, A, S, P. Their perturbative expansion using the MS 
renormalization scheme for the strong coupling reads 

*V) = 4°V) + (^!2) rgV) + (^) 2 rg^foV) + . , (9) 

where the light quark 0(a 2 ) corrections determined in this work will be referred to as r^' MS . As we 
use the pole definition for the heavy mass M throughout this work, only the 0(a 2 ) coefficient of Rq 
is renormalization group /i-dependent. For simplicity we call the quantities = V, A, S, P, decay 
rates from now on in this paper. Because the calculations for the virtual light quark corrections are 
most easily performed by using dispersion integration techniques, the results for all four heavy quark 
currents are derived in the framework of on-shell renormalized QED, where the coupling a is defined 
as the fine structure constant and the wave function renormalization constant of the quark fields is 
fixed by the requirement that the residue of the fermion propagator is one. In a second separate step 
the results are then transferred to the QED MS scheme for the coupling and to QCD. Throughout this 
paper the heavy quark will be denoted as Q with pole mass M and the light quarks as q with mass m. 
It should also be noted that, even in the framework of QED, we will refer to the fermions as quarks 
though they may be leptons as well. For simplicity all quarks are assigned electric charge one. The 
generalization to different charge assignments is straightforward. All QED results are marked with a 
tilde. 

The program of this work is as follows: in Section [2] we describe in some detail the calculation 
of the light quark 0(a 2 ) corrections to massive quark pair production through the vector current in 
on-shell QED. Virtual and real contributions as well as the inclusive sum are discussed separately. 
The results are presented in a general form which is common for all the cases (V, A, S, P). The 
corresponding results for corrections to axial vector induced quark pair production are then given 
in Section ||, whereas in Section || the cases of scalar and pseudoscalar currents are treated on the 
same footing. In Section || we carry out the transition from on-shell QED to the MS scheme and 
to QCD. As mentioned above our results are presented in terms of moments of the inserted vacuum 
polarization from light quarks. In Section |6| this "method of the moments" is discussed in more detail. 
We present the corresponding moments for light scalar and gluonic second order corrections. Section ^ 
contains a discussion of our results in the threshold region and Section |8| our summary. For the sake 
of completeness numerous O{o) formulae for the cases of the different currents (V, A,S, P) are given 
in Appendix [A]. Second order virtual corrections in the case of equal masses m = M are listed in 
Appendix |b|. For the convenience of the reader we also present expansions for the high energy as 
well as for the threshold region in Appendix Phenomenological applications of the results will be 
discussed elsewhere. 
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2 Heavy Quark Production through the Vector Current 



2.1 Virtual Corrections 

The virtual light quark 0(a 2 ) corrections to the massive quark production vertex induced by the 
vector current correspond to the sum of the two particle cuts of the current-current correlator dia- 
grams depicted in Fig. This leads to virtual photon exchange diagrams with an additional vacuum 
polarization insertion of a light quark-antiquark pair. 





Figure 1: Double bubble diagrams where a pair of heavy quarks of mass M is produced through an 
external current. The exchanged photon splits into the light quark-antiquark pair of mass m. 



The basic idea is to write the light quark-antiquark contributions to the photon vacuum polarization 
function U^ u in terms of an one-time subtracted dispersion relation (which accounts for the correct 
subtraction of the subdivergence in the light quark vacuum polarization), 



n^-(mV) 



U q (m 2 ,q 2 



-Ms 



av 2 



Am 2 



A 2 (A 2 



n 9 (mV), 

dA 2 

l 2 - ie) 



Rgq(m 2 ,\ 2 ), 



where 



R qq -(m 2 ,q 2 



1 + 2 



m 



(10) 
(11) 

(12) 



is the absorptive part of the vacuum polarization function of the light quark-antiquark pair, g Thus 
the effective propagator which has to be inserted into the virtual photon line can be written as 



+ ie 

1 

3 



\1<7 



i- 



} 

dA 2 
A 2 



U%(m 2 



+ ie 



pv 



qPq u 



R qq -(m 2 ,\ 2 ) + ... 



(13) 



4 m A 



I 2 — A 2 + ie_ 

The longitudinal terms indicated by the dots are irrelevant due to current conservation. It is evident 



from Eq. (13) that the effective propagator is a convolution of a massive vector boson propagator in 
Feynman gauge with the function R qq -. Inserting expression (13) into the vertex diagrams and carrying 
out the vertex loop integration first, the intended two- loop corrections can be written as a convolution 

2 Via the optical theorem, Eq. ([]), R q q is just the normalized born cross-section for the production of a qq pair in 
e + e~ annihilation into a single photon. 
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integration over the corresponding 0(a) virtual corrections with exchange of a massive vector boson 
of mass A. 

It is customary to parameterize the virtual corrections to the renormalized vector current vertex 
A"j^ in terms of the Dirac (Fi) and the Pauli (F2) form factors 



UQfa) A% v Q (ki) = u Q (k 2 



v Q (h), (14) 



where q = k\ + &2, = § [7^,71/]- uq and vq denote the free heavy quark and antiquark Dirac fields, 
respectively. The perturbative expansions of the form factors read 

F^) = l + ^WV)+(-YVV) + ---, (15) 



7T / \ 7T 



W) = QfP^+Q- FP(^) + .... (16) 

The Dirac form factor is normalized to one at zero momentum transfer q 2 = at all orders in a. To 
0(a) this is completely achieved by the on-shell QED quark field wave function renormalization con- 
stant, which is fixed by the requirement that the residue of the renormalized heavy quark propagator 
is one. No further subtraction is required because the vector current as defined in Eq. (||) is anomalous 
dimension free. Based on the considerations presented above, the light quark corrections to the 0(a 2 ) 
form factors can be written as an one-dimensional integral, 

00 _ 

F%{M 2 ,m 2 ,q 2 ) = \ / F^ (M 2 , q 2 , A 2 ) R q g(m 2 , A 2 ) , i = l,2, (17) 

4 m 2 

where Fy 2 (M 2 , q 2 , A) denotes the first order form factors with virtual exchange of a vector boson 
with mass A. In Eq. (^) we have indicated the mass dependence of the form factors explicitly in 
order to illustrate the mass dependence of the individual factors. It should be noted that the light 
quark corrections to the 0(a 2 ) form factors, Eq. (|T7|), are already renormalized properly and that no 
additional subtraction is needed. In particular, the subdivergence from the virtual light quark loop 
is taken care of properly by using the subtracted dispersion relation in the A 2 integration, whereas 
the overall divergences are subtracted correctly in the 0(a) form factors. The same remark holds 
also true for the corresponding virtual light quark 0(a 2 ) corrections for the axial- vector, scalar and 
pseudoscalar cases presented in the subsequent sections. 
The contributions of the form factors to the decay rate 

iW) = 4°V) + (£) ^V) + ) 2 4V) + ■ • • (18) 



read 



4V) = i/3(3 - /3 2 ) , (19) 



2 



~(l),virt 



r 



V 



(q 2 ) = /3(3-/3 2 ) (ReF 1 (1) (g 2 )+ReF 2 (1) (g 2 ))+/3 3 ReF 2 (1) (g 2 ), (20) 



where 
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is the velocity of the massive quarks in the cm. frame. 



The contribution of F^l to fy 1 is given by 



,(2) 



l%1 



~(2),virt, 
"V,q 



(q 2 ) = (3 - 2 ) ( ReitfJ (q 2 ) + ReF% (q 2 ) ) + 0* ReF% (q 



(2), 



p(2). 



(21) 



The real parts of F^ 2 (M 2 ,q 2 , A 2 ) have already been presented in 14, 18]. For the convenience of the 
reader the full form of the 0(a) form factors with massive vector boson exchange above threshold 
is presented in Appendix A.l| , Eqs. ( |150| ) and ( |151| ). For A 2 — > the expressions for the O(a) form 

are obtained. This results in the well known QED expressions calculated 



factors relevant for r 



(l),virt 
V 



by [|J, where the IR singularity in F\ is regularized by the small fictitious photon mass A, 



F[ L >{M 2 ,q 2 ,\ 2 
1 
2 



1 H j-^ lnp^) In 



1 + /3 2 



1 + 2/3 2 



2/3 



lnp — 1 — i 7T 



M 2 2/3 

1+/? 2 , A 2 ; 1 + (3 2 / 



Li 2 (l-p) + i ln 2 p-3C(2) 



4/3 



1-/3 S 



+ 



1 + 2/3 2 
4/3 



F«(MW,A 2 -0) 



1-/3 2 
4/3 



( Inp 



+ ITT 



(22) 
(23) 



where 



1 + /3 



For arbitrary ratios m/M the evaluation of Eq. (|T^) can be carried out, but is extremely tedious]^] 
However, in the case of interest, m/M — > 0, the task can be simplified enormously by subtracting and 
adding R qq - at its asymptotic value at high energies, R^ = R qq (m 2 , q 2 = oo):[] 



F$(M 2 ,m 2 ,q 2 



dX 2 
A 2 



F^(M 2 ,q 2 ,X 2 ) \R q ,(m 2 ,X 2 



4 m 2 



dA 



+ i R lc I ^F^(M 2 ,q 2 ,X 



A 2 



R q 



1,2. 



4 m 1 



(24) 



Because i?^ reaches its asymptotic value already for values of A 2 much smaller than the heavy scales 
M 2 and q 2 , which govern the O(a) form factors, we can replace F^X(M 2 , q 2 , A 2 ) in the first integral 



on the r.h.s. of (24) by the QED expressions for A 2 — ► 0, Eqs. (^) and (|23"|). The resulting integration 
can then be trivially expressed in terms of the moments Rq and R\ defined by 



Rl 



l 



dx ■ 



ln n x 



R qq (m 2 , 



4 m' 



R q 



n 



0,1,2,... . 



(25) 



The interested reader is referred to [|18|. In the special case m = M the calculation can be simplified considerably 
and leads to the compact results given in Appendix ^ (see also |Q ) . 
4 To our knowledge this trick has been employed the first time in |ig[ |. 
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The evaluation of the second integral on the r.h.s. of (|24|) constitutes the main effort. The final results 
for the light quark 0{a 2 ) corrections to the form factors can be cast into the form (x = q) 



F 



(2) 

1. x 

7 (2) 

2, x 



fi 



(2) 



\n 2T ^--2Rl ln^ + 2^ 

00 q 2 q 2 1 



+ /i 



(l) 



~ Tfl ~ 

R °° 111 ~ R 



+ /l (0) ^, 



r(l) 



where 



fi 



(o) 



fi 



(i) 



fi 



(2) 



R <=° in ~^2 ~ R 



, f (0) nl 



(26) 
(27) 



A virt + -{l + 2(5 2 )B. 



virt 



In 



-l6-) + ^4^ lnp 2^ C(2) + 6 



(28) 



+i 11 { ih 1 2 (1 + p2) ln2 (2 ) " 2 (1 + 2 ^ ln (2) + (1 + ^ c(2) + \ (b + 2 p2) } } 



(1 + (3 2 ) B virt 



1 + 2/? 2 



lnp 



1 + P 2 



+1 7T 



1 + /3 2 



12/3 



1 ln ( l -P 2 ) 

6 V 16 / ' 12/3 ~' 2/3 
2(l + /3 2 ) ln(^)+l + 2/3 2 ]}, 



ln p H + « 7r 

24/3 F 12 



1 + /? 2 
24/3 



C(2) + 



(29) 
(30) 



(o) 



h 
and 



(i) 



(1 — /3 2 ) i?„ ri - 
hip + iir 



3 

1-/3 2 



5 (1 - /3 2 ) 1 1 - (3 2 , s 

' lnp+^-C(2) + in 



12/3 



24/3 

1-/3 2 
12/3 



2/3 



ri-/? 2 




} 


I 6/3 


V2/ 4_ 





A 



virt 



1 + 

6/3 



2 r 



2Li 3 (l-p)+Li 3 (p)+Li 2 (l-p) ln 
1 . , 1 



— hrp+- ln z p ln 2(1- p 



12 



16 



lnp In 2 (2(1 +p) 



2 3 1n(|) + -lnp)C(2)-C(3) 



1 9 
+ — In 2 
12 



1-/3 2 
16 



-Bvirt — . ^ 



1 



Li 2 (l — p) + lnp ( ln(l +p) — - lnp + ln 2 



(31) 
(32) 



(33) 
(34) 



In the real parts of the form factors were already presented in a somewhat more implicit form. 
The light quark-antiquark moments read 



Rio = 


1, 


K = 


ln4--, 

3 ' 


Rl = 


iln 2 4- 
2 


Rl = 


- In 3 4 - 
6 



28 
~9 



28 
9" 



~ ln a 4-| In 2 4+ ( ^ - C(2) ) ln4 - ^ + ~ C(2) + 2C(3) . 



164 5 
~27 + 3 



(35) 
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The reader should note that the entire information on the light quark-antiquark vacuum polarization 
inserted into the photon line is encoded in the three moments R^, Rq and R\. For completeness we 
have also given the moment R%. At this point we would like to mention that the definition of the 
moments of the light quark-antiquark vacuum polarization, as given in Eq. (|2^), strongly relies on the 
dispersion integration approach we used to determine the light quark 0(a 2 ) corrections to the form 
factors. In this approach the occurrence of higher moments R^, n > 1, comes from infrared singular 
terms oc ln n A in the 0(a) form factors for A — > 0|^| If, on the other hand, dimensional regularization 
is employed to regularize ultraviolet as well as infrared divergences, the moments fl35| ) can be uniquely 
identified from the renormalized high energy vacuum polarization function in arbitrary D = 4 — 2e 
dimension. Expanded in terms of small e the vacuum polarization function in the limit q 2 /m 2 — ► oo 
(m being the light quark pole mass and a the fine structure constant) reads (x = q) 



H-fightiq 2 ) — — 3^ 



^0 + R lo ln T— 2 
4m z 



Rf + (R; 



x In— — 



+ R x o 



In- 



\ (R X oo In 



Am 2 



+ R, 



4/i 2 
(in 2 



2 4 m 2 



+ C(2) 



^2 00 4m 2 



Rf) ln 



4 n 2 6 00 4m 2 



+ Ofe 3 



(36) 



where is the common scale parameter introduced in dimensional regularization. Thus, in the frame- 
work of dimensional regularization higher moments i?*, n > 1, in the expression for the second order 
form factors come from (l/e) ra infrared divergences in the 0(a) form factors which cancel the corre- 
sponding e n terms in ^ ight - 



2.2 Real Radiation 

The 0(a 2 ) corrections from real radiation of a light quark-antiquark pair in primary massive quark 
pair production correspond to the sum of all four body cuts of the diagrams shown in Fig. [j] and 
represent the Born rate for primary massive quark production with additional radiation of a light 
quark-antiquark pair. In order to accommodate for the conventions used in the calculation of the 
virtual corrections we will for now remain in the framework of on-shell QED for the calculation of the 
real radiation process. Again we will determine the contribution to the imaginary part of the current- 
current correlator (|]), denoted by fy^ real . This task is accomplished by dividing the four-body phase 
space into three different subprocesses: radiation of a photon with virtuality A 2 off one of the massive 
quarks, propagation of the photon and decay of the virtual photon into the light quark-antiquark 
pair. The first process is described by fy^ real (M 2 , q 2 , A 2 ), the 0(a) contribution to the decay rate 
R v , Eq. @), from real radiation of a photon with virtuality (mass squared) A 2 . The third process is 
proportional to R q q, Eq. ([l2]). Combining both contributions with the photon propagation function, 
summing over all photon polarizations and integrating over all allowed values of A 2 , we arrive at the 

In fact, for secondary production of massive quarks also the moment R\ occurs due to an infrared divergence oc In 2 A 
in the 0(a) form factors describing the primary production of a massless quark-antiquark pair p^| . 
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expression 



M),real 
r V,g 



(M 2 ,m 2 ,q 2 



(y/q 2 -2M) 2 

i r d\ 2 



^r^ real (M 2 ,q 2 ,X 2 )R q ,(m 2 ,X 2 ). 



(37) 



4 m 2 



closely resembles the corresponding formula for the virtual light quark 0(a 2 ) 

2 , q 2 , A 2 ) can be found in Appendix A.l . 

the well known QED result of Q is recovered: 



It is evident that Eq. (|3 
correct 
For A 2 



corrections, Eq. (|17|). An integral representation of r~y^' real (M 



{1) ' real (M 2 ,q 2 ,X 2 ) 



A 2 ^0 



V 



„ j,(2),real 

O Jv 



111 



A 2 



n <-(l),' , eaZ 
JV 



(38) 



where the functions f^' real anc [ jC0> reai are gj ven De l ow ( see Eqs. ( |43| ) and (|44|)). It should be noted 
that the leading infrared behavior for m/M — > of expression (|37|) is determined by the logarithmic 
photon mass singularity in Eq. (|38|). For arbitrary masses m and M the occurrence of three different 
square roots in ([37|) leads to elliptic integrals (regardless of the order of integration) and makes the 
analytic evaluation in terms of polylogarithmic functions (in our opinion) impossible. Even in the 
special case m = M/ the problem of ellipticity persists and no closed result of this four body phase 
space integration is available. This is not surprising as already the massive three particle phase space 
of the first order subprocess cannot be expressed in closed analytic form. However, for any finite ratios 
m 2 /q 2 and M 2 /q 2 the integrations in ( |37|) are well behaved and can be performed in a straightforward 



way numerically [20, 21]. If (at least) one of the quark masses is very light, the problem can be tackled 
analytically: the case of primary production of massless quarks and secondary production of massive 
ones has been completely solved in [12], whereas the analytical result for the opposite mass assignment 
was first presented in In the following we will describe briefly the calculation of the latter and 
give the result in terms of the moments introduced above. 

In contrast to the case of virtual radiation even the calculation of the three body subprocess 
leads to elliptic phase space integrals and could not be carried out explicitly. In the resulting double 
integral (|37|) singularities (in the limit m — > 0) arise from both integrations. In order to deal with this 
complication in a transparent way, we decompose the integration region in soft and hard radiation by 
introducing an infinitesimally small parameter A, m <C A <C M: 



.(2),real 



~{2),real ~(2),real 
1 V, hard "•" ' V, soft 

i r Ai-A/^f 

3 JAM 2 /q 2 



dy + 



(l-2m/^2) 2 



(1-A/y/q 2 ) 



where 



8M 4 , 4M 2 . 
9 4 q 2 



l-y + z)-(l-y + z) 



dy 



2(1 



(i-Vy) 2 dz 



1-y + z 



Am 2 /q 2 



z)y 1 
in — 



— f(y,z)R q q(m 2 ,q 2 z), (39) 

z 



1-y + z 



AM 2 



A 1/2 (y,^) 



l + 



16M 4 , 8M 2 , A 



1 + 



(1-y + z) 2 - (] 



(40) 



with 



l + y 2 + z 2 



2(y + z + yz) . 



(41) 
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In the soft region the energy of the secondary quark pair is smaller than A. The integration area 
is infinitesimally small and hence only the singular part of the integrand contributes. It is therefore 
possible to perform the y-integration first, taking into account only the singular terms. In the region 
of hard radiation, on the other hand, the singularities arising from z — » and 1 — y —* are well 
distinguishable due to our choice A^$> m: 



2A 2m ^_ 

V s V s 



This enables us to separate the singular part easily. By subtracting and adding J-(y, z = 0) we split 
up the integral in two parts: the first part containing [J-(y, z) — J-(y, z = 0)] is infrared-save by 
construction and we can take the limit m — > immediately. The integrand of the second part has a 
much simpler y dependence (through J-(y, z = 0) only) and makes the double integration possible. 

With this decomposition of the phase space in a soft and a hard part the calculation is similar to 
the case of virtual radiation: by subtracting and adding one splits each of the double integrals 
in two parts. In the first part one benefits from the fast convergence of [R q (m 2 ,q 2 ) — R%^\ by taking 
into account only the leading terms of the integrand. The result of the first part is proportional to the 
moments Rq and R\. The second part now becomes tractable as there is one square root less. This 
part is proportional to R%^ and is, similar to the case of the virtual corrections, the most difficult part 
of the calculation.^] 

The result for the rate describing real radiation of a light quark-antiquark in terms of the moments 
is finally given by (x = q) 



-{2), real 
V,x 



r(2),real 
JV 

r(l),real 
JV 



In* 



m 



pi 

2 °° 



rn 



2\n — Rl + 2R 
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TJX 

Kq 



(0),real 
V 



Ri 



+ 



(42) 



with 



r (2), real 
JV 



n(l),real 
J V 



r(0),real 
J V 



3-P 2 
2 

3-P 2 
2 

3-P 2 



a 



real 
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33- 10 P 2 +f3 A 



A r eal ~^ 



15 



48 

6(3 2 
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(3 (39-17/3 2 
24 
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21 



24 

up 2 



U 2 (p) + U 2 (p 2 ) / 

51 + 60/3 + 46 /3 2 



24 

-9 + 24/3 + 24 /3 2 



- Li 2 (p) + 
16 /3 3 - ll/3 4 + 2/3 5 



20 p 3 - lip 4 



+ 



33 - 10/3 2 + /? 4 



+ 



237 



24 

2p 2 - 



24/3 
ln2- 



24 

2 P (-3 + P 2 ) 2 
In p H — In 



(4 - P 2 ) P 2 



96 



51/3 4 P 
Inp H 



ln/3 
-39 + 17P 2 ) 



15 



12 



In 



3 - (l 
22/3 2 + 3/3 4 , (1 - 
24 ln (- 

-P 2 



C(2) 
2p 2 



-P 2 
P 2 



2p 2 



+ 



p (-99 + 37/3 2 
48 



A detailed description of this calculation can be found in 



(43) 
(44) 



Inp 

(45) 
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and 



A 



real 



B 



real 
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real 



1 + /3 2 



{sLi 3 (l -p) + 6Li 3 (p 2 ) + 8Li 3 ( r ^) + 10Li 3 (l -p 2 ) - 13 C(3) 



4 (Li 2 (p)+Li 2 (p 2 )) ln( 



; ' 1 -C(2) [-8 ln2 + 9 lnp + 4 ln/?-6 ln(i— ^) 



2/3 2 



+ 2 lnp ln/3 ( - 2 ln2 + 11 Lap - 16 ln/3-12 In 0— 



4/3 2 

1 - /3 2 

12 In 2 2 lnp + 30 In 2 ln 2 p - 46 ln 3 p- 36 In 2 lnp In ' 



2 



+ 09 hr /, in ( ^-j - 51 lnp In 2 ^—^-) ~ In 3 } . 



1 + /3 2 



Li 2 (p)+Li 2 (p 2 )-2C(2) 



1 /5 

+ - lnpf- lnp + 8 1n2 + 2 1n/3-3 1n(l-/3 2 ; 



2/3 ^ l-/3 2 



()V +^-lnp). 



(46) 



(47) 
(48) 



2.3 Contribution to the Total Rate 

Combining the virtual and real corrections determined in the two previous subsections one arrives at 



_ Jl),virt , ~(l),real 



V 



V 



Z-P 2 B 1-/3 
2 16 



(33 - 39/3 - 17/3 2 + 7/3 3 ) lnp+-/3 (5 - 3/3 ; 



with 



I? 



1 + /3 2 ) (2 1n(l-p) lnp + ln(l+p) lnp + 2 Li 2 (p) + Li 2 (p 2 ) 
-2/3(2 ln(l -p) +ln(l +p)) 

for the 0{a) corrections to Ry, Eq. (|l8[), due to virtual and real photon emission and at 
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A + — [ (l5 - 6/3 2 - /3 4 ) (Li 2 (p) + Li 2 (p 2 )) + 3 (7 - 22/3 2 + 7/3 4 ) Li 2 (p) 
(1-/3) (51 - 45/3 - 27/3 2 + 5/3 3 ) C(2) 
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(1±^1 (-9 + 33 0-9 /3 s -15 s + 4 /J 4 ) hi 2 p 

(33 + 22/? 2 - 7/3 4 ) In 2 -10 (3 - /? 2 ) (l + /? 2 ) ln/3 

1 



(l5 - 22 /3 2 + 3 /? 4 ) ln( 



+ 



2/3 (3-/3 2 ) 
237 + 62 P 2 - 



2 In 

59 /3 4 



2 



ln(l 



i_ 

Af3 2 
P 2 



3 In 



lnp 
1 



/? 2 



4 ln/3 



lnp + 8/9 (3- I3 2 ) (3 ln2-2 ln/3 



6/3 (l7-7^ 



In 



1 -/? 2 ^ /?(75 - 29/? 2 ) 
"2/3 1 



A 



1 + /3 2 



+ 



Li 3 (p)-2Li 3 (l-p)-3Li 3 (p 2 
4 (1 - /3 2 



6 \l+pJ 



5Li,(l 



P 2 ) + yC(3) 



Li 2 (p) In ( - 



+ 2 



Li 2 (p 2 ) In (^) + 2 C (2) ( lnp - In (LJ-) ) 

/? 2 



In 



■1+0 



36 In 2 lnp - 44 ln 2 p + 49 lnp In 



1-/3 5 



+ ln* 



- lnp In (3 (36 In 2 + 21 lnp + 16 In (3 - 22 ln(l - /3 2 )) J 



(52) 



(53) 



for the 0(a 2 ) corrections from virtual and real emission of a light quark- antiquark pair (x = q). It is 
evident that the quadratic logarithms of the ratio m 2 /q 2 in Eqs. (26) and (|4^ ) cancel in the sum (|5l|). 
Further, all terms proportional to R\ add up to zero. As already indicated at the end of Subsection 2A 
this corresponds to the cancellation of the infrared singularity oc ln(A 2 ) in the 0(a) result fy \ when 
the real and virtual contributions are combined. However, one single logarithm of m 2 /q 2 (multiplying 
fy ) remains, rendering the perturbative expansion unreliable if the ratio m 2 /q 2 is very small. This is 
a consequence of the on-shell renormalization scheme, where a is defined as the fine structure constant, 
i.e. at momentum transfer zero. Thus, to determine the correction due to real and virtual radiation of 
massless quarks a mass independent definition of the coupling like the MS scheme is more appropriate. 
The transition to the MS scheme and to QCD will be carried out in Section ||. 



3 Heavy Quark Production through the Axial- Vector Current 



This section is devoted to the presentation of the light quark 0(a 2 ) corrections to the massive quark- 
antiquark axial-vector vertex. As indicated in the introduction only non-singlet contributions are 
considered. We follow the lines of Section |2| and present the results for the virtual and real radiation 
corrections in the framework of on-shell renormalized QED. The transition to the MS scheme and to 
QCD is indicated in Section ||. 

Let us start with the presentation of the virtual corrections to the axial- vector vertex, A^, 
parametrized in terms of the form factors F% and 

7m 75 F 3 {q 2 ) - — (f 75 F 4 (q 2 ' 



UQ{k 2 )K v Q (h) = u Q (k 2 ) 



2M 



(54) 
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where q = k\ + &2- The perturbative expansions of the form factors read 

2 



F*tf) = l+(3^ (1) (. 2 )+g) ifV) + -, 



(55) 



F 4 (q 2 



(56) 



To 0(a) F3 is completely renormalized by the QED quark field wave function renormalization constant. 
In contrast to the Dirac form factor this leads to a non-vanishing value of -F3 for zero momentum 
transfer. Only F3 leads to contributions to the rate R A , 



o 



7T 



f { l\ q 2 )+. 



(57) 



where 



,(0) 



~(l),virt 



P 3 , 

2 (3 3 ReF. 



(i) 



(58) 
(59) 



The contribution of the virtual light quark contributions to f\ is given by 



~{2),virt 
r A,q 



2 /3 3 ReF. 



(2) 



3,9 • 



(60) 



(2) 

In analogy to the calculations for the vector current vertex, Fy 4 q can be written as an integral over 
the 0(a) form factors with exchange of a vector boson with mass A, 



F^(M 2 ,m 2 ,q 2 ) 



1 



dA 2 
A 2 



f[ L \m 2 , q 2 , A 2 ) R q q(m 2 , X 2 ^ 



3,4. 



(61) 



4 m 2 



As in the vector current case Eq. ( pl[) accounts for all necessary subtractions. Explicit formulae for 
.Fg*|(M 2 , q 2 , A 2 ) are presented in Appendix [A.2] for energies above threshold, Eqs. (157) and ( [To" 
Rqq is given in Eq. (|l2|). The well known 0(a) QED expressions are obtained for A 2 — > 0, 



1 



2 + /3 2 



1 + /3 2 , \, A 2 
In p I in — k 



4/5 

F 4 (1) (M 2 ,g 2 ,A 2 



2/3 

lnp — 1 



M 2 



' + ^Li 2 (l-p) + iln 2 p-3C(2) 



2/3 



l TT 



1 + /3 2 A 2 1 + /3 2 
4/3 n M 2+ 4/3 



/1-^n , 2 + /3 2 i 



+ 



4/3 



0) 



(1 - P 2 ) (2 + /? 2 ) 
4/3 



2/3 



(62) 



(63) 



Expressed in terms of the moments introduced in Section ^, Eqs. (|35|), the virtual light quark 0(a 2 ) 
corrections read (x = q) 



5,(2) 

v ( 2 ) 
4,x 



(2) 



r(l) 



2 2 

i?^ ln 2 ^-2i?g ln^ + 2i?f 
2 



(i) 



771 



, f (0) Si 
</ 3 ""-oo ) 



5x 1 m 
K oo m — 2 



+ J4 -n-oo > 



(64) 
(65) 
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where 
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3 
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16 
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V 16 ) 
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i. 



16 

(2 + /3 2 ) ln(^)+i(4 + 3/3 2 ) 

(1 - /? 2 ) (2 + /? 2 ) 



16/3 



4/3 



C(2) + 



2(3 



lnp + 1 



+ ITT 



12/3 



The functions A v i rt and are defined in Eqs. d33| ) and (|3^), respectively. 

The corresponding real radiation corrections can be written as an integral similar to 



(66) 



(67) 



(69) 
(70) 



~(2),real , ,^2 2 2\ 
r A, g (M ,m ,q ) 



(yJq2-2My 



^r^ real (M 2 ,q 2 ,X 2 ) R q ,(m 2 ,X 2 ), 



4 m 2 



(71) 



where f^' rea '(M 2 , q 2 , A 2 ) describes the real radiation of a vector boson with mass A off one of the 
primary heavy quarks. An integral representation of r^' real {M 2 , q 2 , A 2 ) is given in the Appendix 
For A 2 — > the result of [g] is recovered, 

~(l),rea2 -\ 2 ->0 



„ j,(2),real 



In 



f^ 2 
A 2 



q j?(l)ifeai 

-3 Ja 



(72) 



with f^' real anc i j?W> real g[ ven De low in Eqs. (74) and (fT5|), respectively. The complete form of the 
light quark second order real corrections expressed in terms of the moments reads (x = q) 
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/?'"' = P 2 A real + ^f- ( 3 + 2/3 2 + 3/3 4 ) (U 2 (p) + U 2 (p 
1-/3 2 



(0),real _ R 2 A . _,_ 5 ~ P 

96 

(39 + 62/3 2 + 3/3 4 ) Li 2 (p) + ^ (-69 - 37 /3 2 + 160 /3 3 + 33 /3 4 + 9/3 6 ) C(2) 
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48 * ' 



The functions A rea i, B Tea \ and C rea i are denned in Eqs. (fig), (|47|) and fl48|), respectively. 

By adding the corresponding virtual and real contributions one now obtains the 0(a) corrections 
to Ra 

-W = ~(l),viH + ~(l),real (??) 

= /3 2 B-^-^ (21 + 21/3 + 80/3 2 - 16^ 3 + 3/3 4 + 3/? 5 ) lnp+^f (-7 + 10£ 2 + /3 4 
62 V / lb V 

and the 0(a 2 ) corrections due to virtual and real emission of a light quark-antiquark pair (x = i 
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16 



The functions A and B are denned in Eqs. (|53|) and (|50|), respectively. 



4 Heavy Quark Production through the Scalar and Pseudoscalar 
Currents 

This section is devoted to the presentation of the light quark 0(a 2 ) corrections to the massive quark- 
antiquark production rate through the scalar and pseudoscalar currents, Eqs. (||) and (||). As indicated 
in the introduction only non-singlet contributions are considered. We follow the lines of Section ^ 
and present the results for the virtual and real radiation corrections in the framework of on-shell 
renormalized QED. The transition to the MS scheme and to QCD is carried out in Section [|. 

The virtual corrections to the scalar and pseudoscalar vertices are parametrized in terms of the 
form factors S and P, respectively. S and P are normalized to one at the Born level and thus have 
the following perturbative expansions 

S(q 2 ) = 1 +(^ (1) tf)+(|s (2) tf) + ..., (80) 

P(Q 2 ) = l+(£) P^( q 2 )+{^) 2 P^(q 2 ) + .... (81) 

We would like to remind the reader that, as defined in Eq. @, the overall divergences for the scalar 
and pseudoscalar form factors are subtracted by the massive quark wave function renormalization 
constant and the massive quark (pole) mass counterterm. The contributions of the form factors to 
the rates 

= >1> 2 > + (^ 4W 2 > + $!■(«•) + ■ ' ' ( 82 > 

read 



4 0) = /? 3 > rP = P, (83) 



s 

~(l),virt = 2 p3 ReS (l) } ^virt = 2/3Re p(l) (g4) 

for the Born and 0(a) coefficients and 

f {2)Mrt = 2 ^3 Re5 (2) ; f m,mrt = 2(3Re p(2) (g5) 

for the virtual light quark 0(a 2 ) corrections. The integral representation for Sq 2 ^ and pjp is given by 

Af)(M 2 ,mV) = i J ^xW{M\ q \\*)R qq -{ m \\*), X = S,P. (86) 

4 m? 

S ( - 1 \M 2 , q 2 , \ 2 ) and P^\M 2 , q 2 , A 2 ), the 0(a) form factors with exchange of a vector boson with 
mass A, are presented in the Appendix [A. 3; and [A.4| , Eqs. ( |160[) and ( |162D , respectively. For A 2 — > 
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the corresponding QED expressions are recovered, 
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in agreement with || 0. The final expressions for the virtual light quark 0(a 2 ) corrections expressed 
in terms of the moments R^, Rq and R\, Eqs. (|35|), read (x = g) 
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The functions A virt and B virt are denned in Eqs. ( |33| ) and (f^), respectively. 

The integral representation of the 0{a 2 ) corrections to R s l p from the real radiation of a light 
quark pair via photon emission reads 
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4 m 2 



(97) 



where f^cH eal {M 2 , q 2 , A 2 ), the first order corrections due to radiation of a vector boson with mass A, 



S/P 



are given in the Appendices A.3| and A.4 in terms of a one-dimensional integral. For A 2 — > the 0( 
corrections from real photon emission are obtained 
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in agreement with || [7j. The final results for the real light quark 0{a 2 ) corrections are 
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The functions A rea i, B rea i and CV ea ; are denned in Eqs. (fi6|), (|47|) and 
To 0(a) the sums of virtual and real corrections read 
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whereas the light quark (D(a 2 ) contributions to rgj p are (x = q) 
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The functions A and B are denned in Eqs. ( |53| ) and (|50|), respectively. 



5 Transition to MS and QCD 

In Sections || to ^| we have presented results for the light quark second order corrections to the 
absorptive part of the current-current correlators (|l|) and (||) in the framework of on-shell renormalized 
QED. We have seen that in the sum of the virtual and real contributions the quadratic logarithms 
of the ratio m 2 /q 2 cancel as a consequence of the cancellation of the soft photon singularities in the 
sum of virtual and real 0(a) corrections. However, a single logarithm remains which renders the 
perturbative expansion unreliable if the ratio m 2 /q 2 is very small, see Eqs. (pl|), (|T8|) and (108). This 



mass singularity is a consequence of the on-shell renormalization scheme we used so far, where the 
coupling a is the fine structure constant, i.e. defined at zero momentum transfer. However, for very 
small ratios m 2 /q 2 a mass independent definition of the coupling like the MS prescription is more 
appropriate. Because the rates R®, = V, A, S, P, are renormalization group invariant and have 
vanishing anomalous dimensions, the fine structure constant a just has to be expressed in terms of 
the MS coupling in order to perform the transition to the MS scheme.^] 

The 0(a 2 ) relation between the fine structure constant and the running MS QED coupling at the 
scale reads (x = q) 



a Ms(M 2 ) 1 hx , m2 

and leads to 



<* = «ms(^) ( 1 + 2 l n j +0(a 3 m ) (111) 



= ~(o) + ( Ums(v h ~(i) _,_ ^ «ms(m ) \ 2 ~m; 



I) f g) + p y r f)^+ ... (112) 
as the perturbative expansion for the rates R @ , G = V, A, S, P, in the MS scheme, where 
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7 We would like to remind the reader that we keep M to be the massive quark pole mass. 



8 



It should be noted that a^g- is not equivalent to the effective electric charge sometimes used in QED calculations. 
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Rq MS are equal to the corresponding on-shell ones given in Eqs 
between the zero-moments i?Q MS 



It is evident that the logarithmic divergence for m — > has been removed. The moments -R^ MS and 

|35|) , In particular, the equality 
and Rq corresponds to the fact that there is no non-logarithmic 
contribution in the relation between the fine structure constant and the MS coupling, Eq. ( |1 1 1|) . 
The moments -R^ MS and R%^, on the other hand, have to be equal because of renormalization group 
invariance. In fact, the combination a 2 R^ /(3tt) is equal to the contribution of a light quark to the 
0(a 2 ) (3- function of the MS coupling. However, we used different names for the on-shell and the 
MS moments in order to indicate that the MS moments could be directly extracted from the MS 
renormalized vacuum polarization function from a massless quark {x = q) 
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without referring back to the on-shell and dispersion integration definition, Eq. (|25|) f | 

To finally arrive at the corrections to Eq. in the framework of QCD, where the exchanged 
vector bosons in the double bubble diagrams of Fig. [l] are gluons, we have to multiply all previous 
results with the proper SU(3) group theoretical factors, N c = 3, T = 1/2 and Cp = 4/3. We arrive at 
(x = q) 



JO) 
r e 

,(2) 

e,x 



N f (0) 



R X ' MS In 4 



.(1) 



N c C F f^ 



=(1) 



<-) 



e = v,A,s,p. 



(115) 
(116) 



with the moments 
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We would like to remind the reader that all quantities with a tilde refer to QED, whereas those without 
a tilde are QCD quantities. 



6 The Method of the Moments 



Eqs. ([|l]), (f78|), ( |108| ) and ( 116| ) provide a simple and unambiguous method to determine the non- 



singlet light quark second order corrections to the rates i? e and -R , = V,A,S,P: one computes 
the moments of the vacuum polarization function from the light quarks, either through the dispersion 
definition, Eq. (p5|), or directly by the calculation of the vacuum polarization function in the high 
energy limit or for vanishing quark mass, Eqs. (|3^) or ( |114| ), and then inserts the moments into 
Eqs. fl5]]) , (f78|) , (|108| ) and (116) . It is obvious that this "method of the moments" can also be applied to 



determine second order non-singlet corrections from other light particle vacuum polarization insertions 



9 For moments of vacuum polarizations of higher order as well as for higher moments (n > 1) of the one loop vacuum 
polarizations discussed here the situation is more complicated and the on-shell moments are in general different from the 
MS ones. 
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into the photon or gluon line. The only condition is that the absorptive part of the inserted vacuum 
polarization function approaches a constant for large momentum transfer. This is equivalent to the 
occurrence of at most one single power of the logarithm ln(g 2 ) in the renormalized vacuum polarization 
function in D = 4 dimensions. We would like to emphasize that the method of the moments and some 
applications have already been presented before in for the vector current induced rate, R v , using 
exactly the same convention as employed here. However, in view of the new results in this work and 
for the convenience of the reader we will illustrate the method of the moments in the following section 
by applying it to determine the second order non-singlet corrections to the rates from the light scalar 
and the gluonic self energy insertions by identifying the corresponding moments. In pl| the method of 
the moments has even been applied to calculate some third and even fourth order corrections to R v . 
The reader interested in those higher order applications is referred to fll| . The higher order moments 
presented there can be directly transferred to the formulae presented in this work. We would also 
like to mention that in the framework of QED the method of the moments can be used to calculate 
hadronic vacuum polarization corrections by determining the moments of the experimental data to the 
total hadronic cross-section Rhad = cr(e + e~ — > hadrons) / '<J po int via Eq. (|^). This has been exploited 



in [19, EG, and shall not be discussed in this work. 



6.1 Light Scalar Second Order Corrections 

In the framework of on-shell QED the moments of the vacuum polarization of a unit charged scalar- 
antiscalar pair can be easily determined via the absorptive part of the vacuum polarization function 
defined in analogy to Eq. (Ill]), 



RsAm^q') = 1-4^ . (118) 




The moments needed for the 0(a ) corrections to the rates can then be easily determined via rela- 
tion p|, 



R, 



oo 
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4' 



1 2 13 1 

R° = -ln 2 4--ln4+ — --C(2), 

1 8 3 9 4 SV ; ' 

R s 2 = lln 3 4--ln 2 4+f--ic(2)) In4-- + -C(2) + -C(3). (119) 

2 24 3 V 9 4 sv ' J 27 3 v ' 2 SW v ' 

For completeness we have also given the moment R\ and i?f although they are not necessary for the 
light scalar 0(a 2 ) corrections to the total rates. As far as the transition to the MS-scheme and to 
QCD is concerned the situation for light scalars is in complete analogy to the one for light quarks. 
Therefore, all statements and formulae of Section [| are also valid for the scalar case (x = s). 
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6.2 Gluonic Second Order Corrections 



The moments of the gluonic (and ghost) contributions to the 0(a s ) gluon propagator corrections, 
which are needed to calculate the 0{a 2 ) gluonic (and ghost) selfenergy contributions to the rates, as 
illustrated in Fig. [|, can be determined by identifying the corresponding moments in the gluonic (and 




Figure 2: Double bubble diagrams, where the massive quark-antiquark pair is produced by the 
external current and the internal bubble is purely gluonic. The ghost diagrams are not depicted 
explicitly. 



diost) vacuum polarization function as defined in Eq. (|TT 



4 

f - M§-^ 2+ (-H 5 J ln4 J> (120) 

where the gauge parameter £ is defined via the gluon propagator in lowest order 



7 2 + i e 



■<r+c 



Obviously the moments and Tq^' ms , = V, A, S, P, are not gauge invariant. However, for £ = 4 
the combination a 2 R^ AS /(3n) coincides with the complete gluonic and ghost 0{a 2 s ) contributions to 
the QCD /3-function, and thus r@^ MS |^=4 accounts for the leading logarithmic behavior of the sum 
of all gluonic 0(a 2 ) contributions (i.e. proportional to CaCf) in the high energy limit. It is quite 
trivial that such a £ can be found, but it is a remarkable fact that even in the non-relativistic limit, 
i.e. keeping only the leading terms for (3 = y/1 — 4M 2 /q 2 — > 0, rQ^ MS |^ = 4 accounts for all 0{a 2 ) 
corrections to the rates proportional to the color factors CaCf- This issue has already been pointed 
out in [0] for the vector current induced rate and has been used in |1C] for the construction of Pade 



approximants. In Section [7] we will explicitly demonstrate this fact also for the axial-vector, scalar 
and pseudoscalar current induced rates. 



7 Examining the Region for small j3 — close to Threshold 

In the previous sections we have determined the 0(a 2 ) light quark non-singlet corrections to the total 
heavy quark-antiquark production rates induced by the vector, axial-vector, scalar and pseudoscalar 
currents for all energies above the heavy quark-antiquark threshold, q 2 > AM 2 . Using the method 
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of the moments we are further able to calculate the 0(a 2 ) corrections to the rates proportional to 
the color factors CaCf from the (gauge dependent) gluonic and ghost contributions of the 0(a s ) 
corrections to the gluon propagator. 10 Having these two-loop results at one's disposal it is of special 
interest to consider the kinematic regime close to the threshold, where = a/1 — 4M 2 /q' 2 <C 1. Using 
the expansions for — > given in Appendix |C], the 0(a 2 ) light quark non-singlet corrections to the 
rates for all currents read 
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+ 0(/3 3 ) 



where we have displayed an expansion up to next-to-next-to-leading order in j3. The corresponding 
0(a 2 ) contributions from the gluon propagator corrections are straightforward using the gluonic mo- 
ments given in Section 3.2. We would like to emphasize that the expansions (121)-(124) have to 
be considered strictly in the context of fixed order multi-loop perturbation theory which represents 



primarily an expansion in the strong coupling a s . Therefore, the expressions ( |121| )-( |124| ) actually are 
only good approximations in the kinematic regime where a s ^ /3 < 1. It is evident that this kinematic 
regime is rather small, in particular for the cc and and the bb thresholds, where the strong coupling 
becomes quite large. Nevertheless, it is very instructive to examine the expansions because they allow 
for some insights into the long- and short-distance structure of the C(a 2 ) light quark contributions to 
the rates. 

Let us start by comparing the expansions ( |1 2 1| ) — ( p~2^) with the expressions of the 0(a s ) corrections 
in the same limit, 
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The effects of light scalars are straightforward using the method of the moments and will not be discussed further 
in this section. 
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Evidently the expansions ( |121| )-( |l24| ) can be used to carry out an analysis of the respective BLM scales. 
This reveals that the leading and next-to-next-to-leading contributions in the small /3-expansion of the 
0{a s ) contributions to the rates are governed by scales of order 2M(5, the relative momentum of the 
produced heavy quark-antiquark pair, indicating that they are of long-distance origin, i.e. generated 
at non-relativistic momenta. All the next-to-leading contributions for small /3, on the other hand, are 
governed by a hard scale of the order of the heavy quark mass. This clearly shows that the latter 
are of short-distance origin, i.e. generated at relativistic momenta.]^] The most striking feature of the 
expansions ( |121| )-( |i~24| ), however, is that the leading terms in the small /3-expansion for the vector and 
pseudoscalar induced rates, on the one hand, and for the axial-vector and scalar current induced rates, 
on the other, are proportional to each other. Further, the axial-vector and scalar current induced rates 
are suppressed by f3 2 with respect to the vector and pseudoscalar induced rates. For the rest of this 
section we will concentrate only on these dominant contributions for the expansion in j3. In particular, 
we will demonstrate explicitly by using non-relativistic considerations that the latter contributions are 
uniquely determined by the static color potential^ and the orbital angular momentum states in which 
the heavy quark-antiquark pair is produced by the currents. Because the static potential and the 
orbital angular momentum represent all the information needed to completely describe the heavy 
quark-antiquark pair in the non-relativistic limit, the calculations in the following will also illustrate 
the fact that for small j3 the 0(a 2 ) light quark corrections from secondary heavy quark production 
(class(ii)) and the 0{a 2 s ) light quark singlet corrections (class(iii)) are suppressed with respect to 
the non-singlet contributions calculated in this work. Using the same line of reasoning we will also 
determine the leading O^aCfoi 2 ) corrections in the non-relativistic limit and demonstrate that they 
are equal to rQ^ MS |^ = 4, proving the statement given in Section |6.2| . 

For illustration let us first discuss the simpler case of the 0(a s ) corrections to the rates in the non- 
relativistic limit, Eqs. ( |125| )-( |i~28| ). The dominant contributions for small /? can be reproduced by the 
following consideration: the vector and pseudoscalar currents produce the heavy quark-antiquark pair 
in a (relative orbital angular momentum) S-wave state, whereas the axial-vector and scalar currents 
lead to a quark-antiquark pair in a P-wave (see e.g. (^3|). Via the optical theorem the vector and 
pseudoscalar rates are therefore proportional to the absorptive part of the zero-distance S-wave Green 
function, and the axial-vector and scalar rates are proportional to the absorptive part of the zero- 
distance P-wave Green function of the non-relativistic Schrodinger equation 
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G c (x,y) = 5®(x-y), (129) 



where 



a MS (u 2 ) 

V QCD (x) = -C F ^-^ (130) 



11 For the vector and scalar current induced rates these statements can already be found in jlj, |l5| ], 

12 Because the heavy quark-antiquark pair is produced in a color singlet state, we refer only to the color singlet part 
of the QCD potential. 
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is the leading order static QCD potential, with the strong coupling fixed at the scale /i and E = 
\fcp — 2M is the energy relative to the threshold point. Because we are only interested in the regime 
where a s <> (5 <C 1 we can determine G c perturbatively using Rayleigh-Schrodinger time-independent 
perturbation theory (TIPT) starting from the Green function of the free Schrodinger equation, 



M exp(iM /3\x-y\) 
47r \x — y\ 



For the 0(a s ) contributions to the rates it is sufficient to use first order TIPT, 
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The S-wave contribution in Gc at zero distances is just equal to G c (x, 



contribution at zero-distances reads [d xt d y % G c (x, 
over the spatial index i is understood. Taking into account the normalization of the rates as defined 
in Eqs. @ and <S) one finally arrives at 
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which is equivalent to the dominant contributions in Eqs. ( |125| )-( |l2"§| ). Obviously the (3 2 suppression 
of the axial-vector and scalar rates comes from the fact that the latter correspond to heavy quark- 
antiquark pairs in a P-wave state. 

It is now an easy task to determine all C(a 2 ) corrections to the rates in the non-relativistic limit 
coming from light quarks by taking into account the radiative corrections to the static QCD potential 
due to one massless quark species p4j in the configuration space representation pjj 
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Because the corrections to the QCD potential from massless quarks originate completely from the 
Serber-Uehling light quark vacuum polarization [p^| , Eq. ( |114[) , 5V QCD:g can be completely expressed 
in terms of the moments -R^, MS and i?Q MS . Using the same line of reasoning as presented for the 0(a s ) 
contributions to the rates, the 0(a 2 ) light quark corrections in the non-relativistic limit can now be 
determined by calculating the first order (TIPT) correction to the Green function due to 6V QC o,q, 



G^ q (x,y) = - Jd 3 z G^(x,z) 6V QCD>q (z) G^(z,y) 



(138) 
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Extracting the zero-distance S- and P-wave contributions in analogy to Eqs. (|133j )- (|136j ) we arrive at 
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which is in agreement with the dominant contributions in Eqs. ( |121| )-( |l2"4| ). This also illustrates the 
statement that the 0(a 2 ) light quark non-singlet corrections dominate the contributions from classes 
(ii) and (hi) in the non-relativistic limit. Further, Eqs. ( |139| )-( |l42l ) demonstrate explicitly why the 
vector and pseudoscalar rates, and the axial-vector and scalar rates are proportional to each other in 
the non-relativistic limit. 

In the same manner we can now determine all the 0{CaCfch 2 ) contributions to the rate in the 
non-relativistic limit by calculating the first order (TIPT) corrections to the Green function due to 
the gluonic corrections to the QCD potential proportional to CaCf [26], 
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It is a non-trivial fact that the gluonic corrections to the QCD-potential can be expressed entirely in 
terms of the gluonic and gauge dependent moments i?^ MS and i?Q' MS for the special choice of the gauge 
parameter £ = 4 because the CaCf corrections to the QCD potential arise from the combination of 
one-loop propagator corrections, box diagrams and gluonic vertex corrections. As already mentioned 
in Section |6.2| it is quite trivial that such a value of £ can be found for the logarithmic contribution in 
5Vq Cd ^, but it is a remarkable fact that also the non-logarithmic term is accounted for completely by 
i?Q ,MS for the same value of £. Proceeding along the lines of the 0(a 2 ) light quark corrections we can 
now easily determine the 0(CaCfO- 2 ) corrections to the rates in the non-relativistic limit, 
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which proves our statement given in Section 6.2 



To conclude, we would like to note that although the results in Eqs. ( |121|) -( |l24D (and also 
Eqs. (|l39j )- (|l42j )) are only applicable for a s <^ j3 <C 1, they represent very important results for 
predictions of the rates for \(3\ <> a s because they allow for an extraction of the short distance ef- 
fects coming from the light quark dynamics. These short-distance effects are universal for \j3\ <C 1, 
regardless whether (5 is smaller or larger than a s . Therefore, in an effective field theory approach like 
NRQCD D 



, results (121)- (|124| ) are an important input for the so called "matching calculations" 
which allow for a systematic extraction of short-distance contributions. Such a treatment, however, is 
beyond the scope of this paper and shall be carried out elsewhere. 



8 Summary 

The non-singlet 0(a 2 ) corrections from light quarks to the rates for primary production of massive 
quark pairs have been calculated analytically for the production through the vector, axial-vector, 
scalar and pseudoscalar currents. The results have been expressed in terms of moments of the vacuum 
polarization from the light quarks which allows for an immediate determination of C(a 2 ) corrections 
from light colored scalar particles and from the gluonic vacuum polarization insertion. We have 
analyzed the results in the kinematic regime close to the massive quark-antiquark threshold and 
have shown that the dominant contributions in the non-relativistic limit can be reproduced from 
non-relativistic considerations using only the radiative corrections to the QCD potential and the 
information on the orbital angular momentum state in which the massive quark pair is produced. 
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A First Order Corrections from Virtual and Real Radiation of a 
Massive Vector Boson 



In this appendix we present the 0(a) QED corrections to the heavy quark production vertices from 
virtual and real radiation of a vector boson of mass A. Virtual and real radiation corrections are 
displayed separately. The virtual corrections are presented in terms of the functions 
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whereas the real corrections are expressed in terms of an one-dimensional integral representation. The 
parameter t is defined as 

A 2 



M 2 ' 



and 



AM 2 



P 



1-/3 



The real parts of the vector current form factors F^~\ have already been given in [14, 18]. 
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A. 2 Axial- Vector Current Vertex 
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where the functions L and D are given in Eqs. (|154|) and (155). 
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A. 3 Scalar Current Vertex 
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A. 4 Pseudoscalar Current Vertex 
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B Second Order Corrections from Virtual Massive Quarks 

The main focus of this paper is the presentation and discussion of the effects of light quark vacuum 
polarization insertions into the gluon exchange diagrams for massive quark-antiquark production. 
This is not an easy task because for the mass assignment m <C M infrared singularities in the form of 
logarithms of the ratio m 2 /q 2 have to be handled. For the mass constellation m = M, however, it has 
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been observed in [14| that the integrations in Eq. ([l?]) become trivial because the square root connected 
with the born cross-section R q g can be substituted away. The results can then be expressed entirely 
in terms of logarithms. In this appendix we present the results for the axial-vector and pseudoscalar 
second order form factors for the massive quark pair production vertices coming from the the massive 
quark-antiquark vacuum polarization insertion into the gluon line, i.e. for m = M. For completeness 
and convenience of the reader we also review the older results for the real parts of the vector and scalar 
form factors of [|l4|] and p5[ . All the imaginary parts of the form factors are new and have to our 
knowledge never been presented in the literature before. For simplicity all formulae in this appendix 
are presented in the framework of on-shell QED. 
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B.2 Axial- Vector Current Vertex 
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C Expansions 

Below we list the expansions of the quantities Vq' 1 ^ and Sq\ = V, A, S, P, for energies close to the 
heavy quark-antiquark production threshold, (3 = y/l — AM' 2 /q 2 — > 0, as well as in the high energy 
limit, x = M 2 /q 2 -> 0. 
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